The pion parton distribution function in the valence region 
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The parton distribution function of the pion in the valence region is extracted in a next-to-leading 
order analysis from Fermilab E-615 pionic Drell-Yan data. The effects of the parameterization of 
the pion's valence distributions are examined. Modern nucleon parton distributions and nuclear 
corrections were used and possible effects from higher twist contributions were considered in the 
analysis. In the next-to-leading order analysis, the high-a; dependence of the pion structure function 
differs from that of the leading order analysis, but not enough to agree with the expectations of 
pQCD and Dyson-Schwinger calculations. 
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PACS numbers: 14.40.Aq, 25.80.Hp, 12.38.Bx, 12.28.Lg 

The pion has a central role in nucleon and nuclear 
structure. The pion has not only been used to explain 
the long-range nucleon-nucleon interaction, but also to 
explain the flavor asymmetry observed in the quark sea 
in the nucleon. Experimental knowledge of the parton 
structure of the pion arises primarily from pionic Drell- 
Yan scattering |iJ- LJ Li Ll from nucleons in heavy nuclei. 
Recently, the shape of the extracted pion parton distri- 
bution function (PDF) at high- a;, where x is the fraction 
of the pion momentum carried by the interacting quark, 
i.e. Bjorken-a;, has been questioned |5|. 

The anomalously light pion mass is believed to arise 
from dynamical chiral symmetry breaking. Any model 
of the pion must account for its dual role as a quark- 
antiquark system and the Goldstone boson of dynami- 
cal chiral symmetry breaking. Theoretical descriptions 
of pionic parton structure at high-a; disagree. The 
parton model [3, perturbative quantum chromodynam- 
ics (pQCD) 0, 11] and some non-perturbative calcula- 
tions such as Dyson-Schwinger Equation (DSE) mod- 
els 0, 0, Lul LLH El indicate that the high-a; behavior 
should be near (1 — x) a , with a > 2 . In contrast, rela- 
tivistic constituent quark models |13L Il4| , Nambu- Jona- 
Lasinio models with translationaly invariant re gula riza- 
tion d, EE E3, , Drell-Yan- West relation Haljj and 
even duality arguments [2l| favor a linear high-x depen- 
dence of (1 — x). Instanton-based models appear to lie 
in between these two [22J . Lattice calculations yield only 
the moments of the distributions and not the PDFs them- 
selves miii. 

While the PDFs for the nucleon are now well deter- 
mined by global analyses of a wide range of precise data 
(see e.g. |25l I2H |27| ) the pion PDFs are very poorly 
known. Presently there are only two experimental tech- 
niques to access quark distributions in the pion: the 
deep inelastic scattering from the virtual pion cloud of 
the proton for which data is available in the low 
region (3 x 1(T 4 < x < 0.01) [M E3, and the Drell- 
Yan mechanism which provides data in the valence re- 
gion (0.2 < x < 0.99). Unfortunately there is no overlap 
between these two experimental techniques. 

A previous leading order (LO) analysis of pionic Drell- 
Yan data by J. Conway et al. (Fermilab E615) [l| sug- 



gested a high-a; dependence of (1 — a;) 1 ' 26 . The data were 
at relatively high momentum transfer, Q 2 > 16 (GeV) 2 , 
where pQCD should be a valid approach. Previous 
analyses in both LO Q, and next-to-leading order 
(NLO) 0, 113 have adopted a simple functional form 
for the pionic valence quark distributions or made other 
assumptions about the parton distributions. Since the 
original analysis of the Drell-Yan data, it has been sug- 
gested Q that the simple functional form assigned to 
the valence quarks in the pion and that using only a LO 
analysis might introduce a bias at high-a;. 

The purpose of this work is to re-analyze the pionic 
Drell-Yan cross section data in order to study the high- 
x behavior of the pion PDFs and to see if such a bias 
exists. Cross section data from Fermilab E615 (Con- 
way et al. 0) were fit to determine the form of the 
high-a; pion parton distribution. A NLO analysis using 
the well understood PDFs for the proton (MRST98 (26|, 
GRV98 il, and CTEQ5M @, including nuclear cor- 
rections [23) and allowing for higher twist effects was 
performed. 

The LO Drell-Yan cross section for a pion interacting 
with a nucleon is 
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where the sum is over quark flavor, q^^) is the PDF 
for quark flavor q in the pion (nucleon); e is the charge 
of the quark, M 1 is the mass of the virtual photon and 
^(N) is the momentum fraction (Bjorkcn-a;) of the inter- 
acting quark in the pion (nucleon). (Where not needed 
for clarity, the subscript 'V will be dropped.) The cross 
section per nucleon on an atom with atomic number Z 
and atomic mass A is 
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Charge symmetry {e.g. M pro t. 



etc.) was then 



used to express the cross section in terms of just the pion 
and proton's PDFs. 
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The pion's valence (<j£), sea(q£) and gluonic(p 7r ) parton 
distributions were parameterized as 
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respectively. The valence parameterization, xq^, follows 
that suggested by Hecht et al. with the addition of 
the term proportional to ar/m?„ that allows for the pos- 
sibility of higher-twist effects as suggested by Berger and 
Brodsky [34J and used by Conway et al. 0,12. The co- 
efficient of this term, 7, was restricted to only positive 
values. This parameterization reduces to the "minimal 
parameterization" used in earlier works Q, Q if e and 
v are both fixed at 0. If (3 = 1, the pion valence PDF 
has a linear high- a; behavior and as (3 increases, there 
is more curvature at high- a;. The pion's sea and gluon 
distributions are identical to those used in earlier works. 
We assume the pion's valence and sea distributions are 
SU(3) flavor symmetric: 



xq^(x) = xu^(x) = xd^(x) and 

xql(x) = xul(x) = xu%{x) 

= xa%(x) = xd B 7r (x) 

= X K( X ) = xsl{x). 



(6) 
(7) 



The sea and gluonic distributions are better determined 
from other data. The shape of the gluon distribution is 
deduced from CERN WA70 prompt photon data [j^]. A 
fit to these data is it best described by 77 = 2.1 and = 
J Q 1 xg^(x)dx = 0.47 [p. The sea distributions can be 
determined by comparing tt + and ir~ induced Drell-Yan 
scattering as was done by the CERN NA3 collaboration 
who found S = 8.4 These parameters were fixed at 
the values given above and were also used by Conway et 
al. 0. 

Sum rules constrain the normalization coefficients 
and A\. Specifically, the total number of valence u (and 
d) quarks in the tt~ is constrained to be unity by 



/ ql (x^dx^ = 1. 
Jo 



(8) 



Momentum conservation within the pion is enforced by 
requiring 



xql (av) dxv +6 xq% (x n ) dx„ +G % = 1. (9) 



Previous experiments have shown that the basic, LO 
Drell-Yan cross-section formula fails to explain the mag- 
nitude of the observed cross-section, cr cxp , by a factor 
of nearly two. This is traditionally accounted for by 
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FIG. 1: (Color online) The ratio of NLO to LO pionic Drell- 
Yan cross sections (_K"-factor) calculated using CTEQ5M pro- 
ton PDF and the pion PDF from this work. One point is 
shown for each (xf, s/t) bin projected onto the NLO 
terms are clearly important asu 1. 



multiplying the LO cross section, er LO by a "if'-factor, 
i.e. a cxp = Ka LO , which is completely independent 
of kinematics. This difference is attributed to NLO 
terms and in proton- nucleon Drell-Yan, K' w 1, where 
CT c Xp = K ' a ^o ;Ui Comparisons of LO and NLO cross 
section calculations show that the "if" -factor is not com- 
pletely independent of kinematics, in particular in the 
high- a; region, as shown in Fig. ^ Note that K and K' 
are indistinguishable from normalization errors. In the 
NLO fits, it was assumed that the NLO calculation ac- 
counted for the entire observed cross section, and so K' 
reflected the normalization uncertainty. Thus, deviation 
of K' from one contributed to the x 2 of the NLO fits, 
while K was allowed to vary freely in the LO fits. 

The differential cross-section measured by Conway et 
al. were the basis for this work. Events arising from 
J/ip and T decays were removed by requiring that 4.05 < 
M 1 < 8.53 GeV. (Relaxing this condition to include 
bins with M 1 > 11.06 GeV had no noticeable effect on 
the results.) Only bins with xp > (Feynman-a;) were 
included in the fit. Furthermore to minimize the effects of 
shadowing in copper, only bins in which x^ > 0.06 were 
used. Finally, if the center of the bin was not within 
the defined acceptance that bin was not considered by 
the fit. It is important to note that this final restriction 
only removed six bins; however, three of these had x v > 
0.97. These same restrictions were used in the original 
fit of these data 0,0. With these restrictions, data from 
76 bins of the 168 tabulated in Ref. 1] remained. The 
averaged mass for the data included in the fit was 5.2 
GeV. 

The data were fit in both NLO and LO. In both cases, 
no QCD evolution was performed so that the resulting 
pion PDFs are at the average mass scale of the data, 
(M 7 ) 2 = (5.2 GeV) 2 . The LO fits were used to com- 
pare with the fit of Conway et al. and so that the 
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TABLE I: Tabulated below are the results from the leading order (LO) and next-to-leading order (NLO) fits at a mass scale of 
{M-y} 2 = (5.2 GeVY^. The original LO fit of Conway et al. [jj is tabulated in the first column. The second column is a LO fit using 
the Conway et al. [IJ parameterization of the proton. Succeeding three columns are based on the CTEQ5L parameterization of 
the proton showing the effects of nuclear corrections and the extended pion PDF parameterization. Following the LO results are 
the results of the NLO fits. The sixth column gives the results for the full parameterization while the seventh column shows the 
"minimal" parameterization results (preferred fit). The next two columns give the results for the "minimal" parameterization 
using the MRST and GRV proton parameterizations. In the final column we see the effect of eliminating the higher twist 
contributions. In all cases, 8 and were fixed at 8 = 8.40 and = 0.47. 







Leading Order Fit Results 




Next-to-Leading Order Fit Results 


Prot. PDF 


Fit from 


Conway 


CTEQ5L 


CTEQ5L 


CTEQ5L 


CTEQ5M 


CTEQ5M 


MRST98 


GRV98 


CTEQ5M 


Nucl. Corr. 


Conway [1] 






EKS98 


EKS98 


EKS98 


EKS98 


EKS98 


EKS98 


EKS98 


a 


0.60±0.03 


0.65±0.07 


0.67±0.07 


0.67±0.07 


0.43±0.17 


0.43±0.05 


0.70±0.06 


0.69±0.06 


0.69±0.06 


0.66±0.04 


P 


1.26±0.04 


1.33±0.08 


1.35±0.08 


1.36±0.08 


1.38±0.26 


1.60±0.08 


1.54±0.08 


1.54±0.08 


1.53±0.08 


1.46±0.04 


e 


0.00 


0.00 


0.00 


0.00 


-3.39±7.02 


-3.30±1.90 


0.00 


0.00 


0.00 


0.00 


V 


0.00 


0.00 


0.00 


0.00 


-0.76±5.52 


-0.01±1.35 


0.00 


0.00 


0.00 


0.00 


7 


0.83±0.26 


0.78±0.65 


0.75±0.63 


0.77±0.62 


3.0±2.8 a 


3.0±2.0 a 


0.60±0.34 0.57±0.34 0.57±0.34 


0.00 


K (K') 


1.75±0.13 


1.57±0.11 


1.53±0.10 


1.49±0.10 


1.56±0.24 


1.02±0.05 


0.97±0.06 


0.98±0.06 


0.98±0.06 


1.01±0.05 


X 2 /Ndf 


359/329 


67.6/72 


67.1/72 


69.3/72 


69.1/70 


72.8/70 


73.1/72 


72.0/72 


70.9/72 


74.4/73 



"Converged at upper limit of parameter. 



effects of the NLO terms in the cross section could be 
clearly understood. The results shown in Tab. |U All 
of the LO fits have slightly more curvature at high- a; 
than the original fit of these data 0. This includes the 
one that used the Conway parameterization of the pro- 
ton. The original work by Conway et al. was based on 
a finer binning of the data, but the largest differences 
between the two fits are in the if-factor. In this con- 
text, it should be noted that a direct calculation based 
on Conway et aVs parameterizations of the proton and 
pion reveal that K = (Jnlo/vlo ~ 1.5 in close agree- 
ment with the present fit and a kinematic dependence 
similar to that in Fig. 2] Fixing K = 1.75, more closely 
reproduces Conway et al.'s fit with a = 0.56 ± 0.01 and 
f3 = 1.24 ± 0.05. The use of modern parton distribu- 
tions and the inclusion of nuclear corrections have only a 
slight effect on the high-x PDF. In LO when e and v were 
allowed to freely vary there were significant correlations 
between fit parameters with e becoming large and nega- 
tive while 7 became large and positive for trivial gains in 
X 2 - A bound of < 7 < 3.0 was implemented in the fit. 
The shallowness of the x 2 hyper-surface is evident in the 
large uncertainties in these parameters. 

These fits were repeated in NLO using the CTEQ5M 
proton parameterization. As expected from the kine- 
matic dependence of the if-factor (see Fig. there is 
even more curvature at high-x in both the "minimal" 
and full parameterization with (3 = 1.54 and j3 = 1.60 
respectively. As in the LO case, the fit to the full param- 
eterization was able to offset increases in 7 (the higher 
twist term) with changes in e and v for extremely slight 
gains in \ 2 - Based on this, and the relatively large un- 
certainties in these three parameters, it was reasonable 
to remove either 7 from the fit or fix e = and v = 0. 
Fixing 7 at values between 0.0 and 0.8 yields fits which 



have significantly larger uncertainties in e and v. Alterna- 
tively, 7 has a direct physical interpretation, while e and 
v dilute the interpretation of (3 and are merely present 
to allow for a better representation of the data. While 
not considered in this work, the angular distributions ob- 
served by Conway et al. and by S. Falciano et al. (CERN 
NA10) |37| are better described with the inclusion of the 
higher twist term. Using the "minimal" parameteriza- 
tion (fixing e = and v = 0) and allowing 7 to vary fits 
the data as well as allowing all three parameters to vary. 
In fact, the distributions of residuals as a function of x 
are nearly identical to the full fit. Based on these con- 
siderations, the "minimal" parameterization was adopted 
as the preferred parameterization. This fit, compared to 
data is shown in Fig.|2and the resulting pion valence par- 
ton distribution in Fig. To facilitate the comparison 
of this fit with lattice results, Tab ITT1 gives the moments, 

x n q^(x)dx, of the pion's valence distributions. 

To determine if there was any bias in the fit, it is useful 
to examine the residual or pull of each data point defined 
by R = (cr mcas - - cr calc ) /<5 meas - where cr meas -( calc -) is the 
measured (calculated) cross section and (5 data is the un- 
certainty of the measured cross section. The distribution 
of the residuals as a function of x reveals no bias, shown 
in Fig. 2] The overall distribution is well described by a 
Gaussian distribution with a standard deviation consis- 
tent with unity and mean with zero as shown in Fig. 

For comparison of the influence of the proton PDF on 
the fit, additional NLO fits were performed using the 
MRST [|(| and GRV98 H3 parameterizations. These 
fits produced results similar to those obtained using the 
CTEQ parameterization, as shown in Tab. [i] 

The contribution of the pion's strange sea was also 
investigated. The strength of the pion's strange sea was 
varied from being equal to the up and down sea [SU(3) 
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FIG. 2: (Color online) Measured values of da/dxFd^/r (points) reported by Conway et al. PJ along with the NLO fit to the 
data (curves). 
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(Color 



online) (a) The residuals, R = 
' s ', for the preferred fit as a function of 
x. No bias is apparent as a function of x. (b) The resid- 
ual distribution fit to a Gaussian with (R) = 0.14 ± 0.15 and 
a = 1.09 ±0.14. 



TABLE II: The moments of the pion valence PDF, evaluated 
at a mass scale of (M 7 ) 2 = (5.2 GeV) 2 , are tabulated below. 
The th moment was a constraint of the fit-see Eq. |H] 
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I n 

/ x q, 
Jci 


/ T (x)dx 



1 0.217±0.011 

2 0.087±0.005 

3 0.045±0.003 



FIG. 3: (Color online) The valence pion PDF, xq^(x), for the 
preferred fit (solid) as well as those from Conway, et al. Q 
(dashed), SMRS 31] (dotted) and GRS 27] (dot-dashed) at 
a mass scale of (M 7 ) 2 = (5.2 GeV) 2 are shown. The present 
work shows evidence for greater curvature at high- a;. The 
upper plot shows the entire x range of the data, while the 
lower plot emphasized the high- a; region. For comparison, the 
small higher twist (7) term has been removed. 



flavor symmetric limit-see Eq.0 to no strange sea. Since 
most of the data was in the pion's valence region, it is 
not surprising that the strange sea of the pion had little 
effect on the valence parameterization. 

Several earlier works have also considered the data fit 
by the present work, as shown in Fig. [21 The fit by Con- 
way et al. was in LO and showed less curvature at high-x, 



5 



but Conway even speculated that a NLO fit might yield 
different results @. The fit by SMRS H3 looked at this 
data as well as other pionic data. While this fit was in 
NLO, it was based on what is now an old description of 
the proton's PDFs with a linear, ad hoc nuclear correc- 
tion. In addition, some of the high-r data was excluded 
from the fit. GRS H3 studied the pionic PDF in NLO 
with the less restrictive parameterization, but this work 
also invoked "constituent quark independence" which di- 
rectly relates the pion and proton distributions |27| . in- 
troducing a different possible bias. 

We have presented a NLO analysis of Drcll-Yan data 
from Fermilab E-615. When analyzed in NLO, the pi- 
ons valence PDF clearly have more curvature as x — ► 1 
than previous leading order fits. This is in worse agree- 
ment than previous analysis with the linear dependence 



expected by NJL models [H E IH 13 , duality argu- 
ment H2 and the Drell-Yan-West relation 0, Ojij. It 
is not as much as is expected from Dyson- Schwinger 
Equation based models of the pion 0,0, 0, ITTl fl^ | 
or pQCD [71,13, however. Finally, while E615's resolu- 
tion in x is not given, the mass resolution for the J /ip is 
o~m = 0.16 GeV Q. This could translate into substan- 
tial uncertainty in x as x — > 1, so additional data from 
a more precise experiment with better resolution would 
certainly be welcome. 
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